MAXIMUM PRINCIPLES FOR PI-CONFORMING FINITE 
ELEMENT APPROXIMATIONS OF QUASI-LINEAR SECOND 
ORDER ELLIPTIC EQUATIONS 
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Abstract. This paper derives some discrete maximum principles for Pl-conforming finite el- 
ement approximations for quasi-linear second order elliptic equations. The results are extensions 
of the classical maximum principles in the theory of partial differential equations to finite element 
methods. The mathematical tools are based on the variational approach that was commonly used 
in the classical PDE theory. The discrete maximum principles are established by assuming a prop- 
erty on the discrete variational form that is of global nature. In particular, the assumption on the 
variational form is verified when the finite element partition satisfies some angle conditions. For the 
general quasi-linear elliptic equation, these angle conditions indicate that each triangle or tetrahe- 
dron needs to be 0(h a )-a.c\ite in the sense that each angle ctij (for triangle) or interior dihedral angle 
ctij (for tetrahedron) must satisfy onj < it/2 — ■yh a for some a > and 7 > 0. For the Poisson prob- 
lem where the differential operator is given by Laplacian, the angle requirement is the same as the 
existing ones: either all the triangles are non-obtuse or each interior edge is non-negative. It should 
be pointed out that the analytical tools used in this paper are based on the powerful De Giorgi's 
iterative method that has played important roles in the theory of partial differential equations. The 
mathematical analysis itself is of independent interest in the finite element analysis. 
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1. Introduction. In this paper we are concerned with maximum principles for 
PI conforming finite element solutions for quasi-linear second order elliptic equations. 
The continuous problem seeks an unknown function with appropriate regularity such 
that 

(1.1) - V ■ (a(i, 11, Vu)Vm) + b(i, «, Vm) ■ V11 + c(i, «)w = /(a;), in f2, 

where 51 is a polygonal or polyhedral domain in R d (d = 2,3), a = a(x,u,Vu) is 
a scalar function, b = (bi(x 7 u,\7u))dxi is a vector-valued function, c = c(x, u) is a 
scalar function on tt, and Vu denotes the gradient of the function u — u(x). We shall 
assume that the differential operator is strictly elliptic in ft; that is, there exists a 
positive number A > such that 

(1.2) a(x,T],p) > A, Vx 6 n,T] G R,p G R d . 

We also assume that the differential operator has bounded coefficients; that is for 
some constants A and v>0we have 

(1.3) \a(x,r],p)\ < A, A~ 2 ^ \h(x, r/,p)\ 2 + A~ 2 |c(ir, r/)\ 2 < v 2 , 
for all x G ft, 77 G M, and p G R d . 
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Introduce the following form 

(1.4) Q.(w;u,v) := / {aVti • Vw + b • (Vm)w + cuv} dx, 

Jn 

where a = a(x, w, Vw), b = h(x, w, Vw), and c = c(x, iw). Let the function / in (|1. II) 
be locally integrable in f2. Then a weakly differentiable function u is called a weak 
solution of (jl.ip in fl if 

(1.5) Q(«;«,w)=F(«), VveCg(fi), 

where F(u) = / n fvdx. For simplicity, we shall consider solutions of (jl.ip with a 
non-homogeneous Dirichlet boundary condition 

(1.6) u — g, on 90, 

where g £ H^(dfl) is a function defined on the boundary of f2. Here i? 1 (SI) is 
the Sobolev space consisting of functions which, together with its gradient, is square 
square integrable over ft. H^(dft) is the trace of H 1 (ft) on the boundary of ft. The 
corresponding weak form seeks u £ H (ft) such that u = g on dfl and 

(1.7) Q(u;u,v) = F(v), Vv € H^ft). 

The usual maximum principle for the solution of (|1.7j) (e.g., see [9]) asserts that 
if c(x, rf) > and f(x) < for all x G f2 and 77 € M, then 

(1.8) supu(a;) < sup g+(x), 

xeQ xedfi 

where g+(x) = max(g(x), 0) is the non-negative part of the boundary data. Moreover, 
if c = 0, then one has 

(1.9) supu(a;) < sup g(x). 

xeti xedfi 

For general non- homogeneous equation (jl.ip . by using the powerful De Giorgi's it- 
erative technique [6] one can derive the following maximum principle (see |20j for 
details). 

Theorem 1.1. Let u € H 1 ^) be a weak solution of U.l\) and il.6)) arising from 
the formula \1. 7| ). Let p > 2 be any real number such that p < +00 for d — 2 and 
p < ford > 2. Assume that f 6 L (p-ih^-i) (ft) with a real number r, 1 < r < p—1. 
Assume that the coefficient functions and the solution satisfy c — • b > for any 
x G ft. Then, there exists a constant C = C(fi) such that 

(1.10) supu(x)<k + C\\f\\ LT _^ TJ 
where 



fc 



sup g+(x), if c > 0, 

xGdQ 

sup g(x), if c = 0. 



Moreover, the dependence of C = C(ft) is given by 

C(ft) = C2t^\flfl£ 
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Here and in what follows of this paper, C denotes a generic dimensionless constant. 

The goal of this paper is to establish an analogy of the maximum principles (|1.8p , 
p. 91) . and f| 1 . 10[) for Pl-conforming finite element approximations of (ll.7p . We will 
establish similar maximum principles for such finite element approximations with an 
assumption on the form Q[w; u, v) (see (|3.6p for details) that can be verified through 
some geometric conditions imposed on the corresponding finite clement partition. As 
an example, we shall explore some geometric conditions that apply to the angles 
of each element, as was commonly done in existing results on discrete maximum 
principles (DMP) (see for example, [5] and [18]). For the general quasi-linear elliptic 
equation (jl.ip . the triangles or tetrahedron need to be 0[h a )-&cute in the sense that 
each angle (for triangular case) or interior dihedral angle (for tetrahedral case) must 
satisfy Oij < tt/2 — ^h a for some a > and 7 > 0. For the Poisson problem where 
the differential operator is given by Laplacian, the angle requirement is the same 
as the existing ones: either all the triangles are non-obtuse or each interior edge is 
non- negative as defined in |S]. 

The research on discrete maximum principles for finite element solutions can 
be dated back to the seventies of the last century. In [5], a linear second order 
elliptic equation was considered, and a discrete maximum principle was established 
for continuous piecewise linear finite element approximations if all angles in the finite 
element triangulation are not greater than tt/2 (the so-called non-obtuse condition). 
In [TS], it was noted (see page 78) that the discrete maximum principle holds true 
for continuous piecewise linear finite element approximations for the Poisson problem 
under the following weaker condition: for every pair [0.1,0.2) of angles opposite a 
common edge of some given pair of adjacent triangles of the triangulation one has 
o\ + 02 < 7T. In 16 , it was shown that the discrete maximum principle may hold 
true in some cases if both angles in such a pair are greater than tt/2. In [3], the case 
of rectangular meshes and bilinear finite element approximations was considered for 
second order linear elliptic equations with Dirichlet boundary conditions. The notion 
of non-narrow rectangular element was introduced as a sufficient geometric condition 
for a discrete maximum principle to hold. In |14j . a 3D nonlinear elliptic problem 
with Dirichlet boundary condition was considered and the effect of quadrature rules 
was taken into account. A corresponding discrete maximum principle was derived 
under the condition of non-obtuseness for the underlying tetrahedral meshes. It was 
further shown that the DMP may also hold true for continuous piecewise linear finite 
element approximations for elliptic problems under various weaker conditions on the 
simplicial meshes used. The acuteness assumption has been weakened in [13| and 
[16] . In particular, in certain situations, obtuse interior angles in the simplices of 
the meshes are acceptable. In [llj . quasi-linear elliptic equation of second order 
in divergent form was considered, and corresponding DMPs were derived for mixed 
(Robin- type) boundary conditions. In [T7], a weaker discrete maximum principle is 
shown to hold under quite general conditions on the mesh (quasi-uniformity) and 
arbitrary degree polynomials, namely 

Halloo, n < C\\u h \\ 

where C > is independent of the meshsize h. In [8], positivity for discrete Green's 
function was investigated for Poisson equations. The authors addressed the question 
of whether the discrete Green's function is positive for triangular meshes allowing 
sufficiently good approximation of H 1 functions. They gave examples which show 
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that in general the answer is negative. The authors also extended the number of cases 
where it is known to be positive. 

The contributions of this paper are as follows: (1) the DMP result with general 
non- homogeneous quasi-linear elliptic PDE ([l.ip is new (see Theorem 13. 2|) : (2) the 
DMP result, as summarized in Theorem 13. 3[ is new with the inclusion of the first 
order term h(x, u, Vu) ■ Vu in the PDE; and (3) the mathematical tools for deriving 
DMPs are new in the finite element analysis. Our analytical tools are based on a 
variational approach which are extensions of similar tools that were used to derive 
maximum principles in pure theory of partial differential equations. We envision that 
the new analytical tool shall have applications to a much wider class of problems than 
the existing approach based on the inversion of M-matrices in the DMP analysis. 
In particular, we shall report some DMPs for Pl-nonconforming finite elements and 
mixed finite element approximations for (jl.ip and (|I.6[) in a forthcoming paper. 

The paper is organized as follows. In Section 2, we shall review the finite element 
method for (jl.lj) and (|1.6[) based on the form (|2.3[) . In Section 3, we shall derive two 
discrete maximum principles (DMP) for Pl-conforming finite element approximations 
under an assumption to be verified in forthcoming sections. In Section 4, we discuss 
the relation of shape functions with angles and interior dihedral angles for each element 
(triangular or tetrahedral). Finally in Section 5, we shall verify the assumption under 
which the DMPs were derived in earlier sections by requiring some angle conditions 
for the underlying finite element partition. 

2. Galerkin Finite Element Methods. In the standard Galerkin method 
(e.g., see [H H]), the trial space iP(fi) and the test space Hq(H) in (|1.7[) are each 
replaced by properly defined subspaces of finite dimensions. The resulting solution 
in the subspace/subset is called a Galerkin approximation. Galerkin finite element 
methods are particular examples of the Galerkin method in which the approximating 
functions (both trial and test) are given as continuous piecewise polynomials over a 
prescribed finite element partition for the domain, denoted by Th- 

We consider only Galerkin finite element approximations arising from continuous 
piecewise linear finite element functions - known as PI conforming finite element 
methods. To this end, let Th be a finite element partition of the domain Q consisting 
of triangles (d = 2) or tetrahedra (d — 3). Assume that the partition Th is shape 
regular so that the routine inverse inequality in the finite element analysis holds true 
(see [1]). Denote by h = max^g^ the meshsize of Th with being the diameter 
of T. For each T E Th, denote by Pj{T) the set of polynomials on T with degree no 
more than j. The PI conforming finite element space is given by 

(2.1) S h := {v : v G H\n), v\ T G Pi(T),VT G T h ) ■ 

Denote by S® the subspace of Sh with vanishing boundary values on dQ; i.e., 

(2.2) S° h :={veS h , v\ dn =0}. 

The corresponding Galerkin method seeks Uh G Sh such that Uh — Ih9 on <90 and 

(2.3) Q(u h ;u h ,v) = F(v), Vu G 5°, 

where Ihg is an appropriately defined interpolation of the Dirichlet boundary condition 
(|1.6[) into continuous piecewise linear functions on dtt. For example, the standard 



5 



nodal point interpolation would be acceptable if the boundary data u = g is sufficiently 
regular. 

Let v € Sh be any finite element function and k be any real number. We shall 
decompose v — k into two components 

(2.4) v - k = (v - k) + + (v - fc)_, 

where (w — fc)+ is a finite element function in SVi taken as the non-negative part of 
v — fc at the nodal points of the finite element partition Th', i.e., (v — fc)+ is defined as 
a function in 5Vi such that at each nodal point A, 



(v-k) + (A) 



v(A)-k, ]£v(A)>k, 
0, otherwise. 



Likewise, the function (v — fc)_ := (v — k) — (v — fc) + is the non-positive part of v — k 
at the nodal points of Th- 

Lemma 2.1. £e£ u £ Sh be any finite element function. Let k be any real number 
such that k > if c — c(x, r) > anrf k arbitrary if c = 0. Then, we have 

(2.5) «, (« - > Q(v; (v - (w - k)+) + 0(«; (« - *)_, (« - *)+). 



Proof. Observe that u) is bilinear in terms of u and v. Thus, 

Q(v; v, (v - k) + ) = Q(v; v-k,(v- k)+) + Q(v; k, (v - k)+) 
= Q{v; v-k,{v- k)+) + k(c, (v - k)+). 

Here we have used the fact that Q(v ; k, (v — fc)+) = k(c, (v — k) + ). If c > and k > 0, 
then we obtain 

(2.6) Q{v; v, (v - k)+) > Q.(v; v — k,(v — k) + ). 

In the case of c = 0, (|2 . 6[) clearly holds true for any real number k and the inequality 
can be replaced by equality. It follows from (|2.6j) and the decomposition (|2.4|) that 
(|2.5|) holds true. This completes the proof of the lemma. □ 

For convenience of analysis, we shall need a discrete equivalence for the usual L p 
norm ||u||lp in the finite element space Sh- To this end, let v be any finite element 
function in Sh- Denote by {v} the vector 

{v}={v(A 1 ),...,v{A j ),...,v(A N )), 

where {^4.j}j=i,— .n is the set of nodal points of the finite element partition Th- Denote 
by Slj the macro element associated with the nodal point Aj (i.e., Vtj is the union of 
elements that share Aj as a vertex point). It is not hard to show that there exist 
constants Co and C\ such that 

N N 

(2.7) C„5>(^)H%I < < C^HAj^m- 

3 = 1 3 = 1 



G 



For completeness, let us outline a proof for the left inequality. For any x £ Clj, we 
have 

v(Aj) = v(x) + (Aj - x) ■ Vv. 

Thus, 

\v(A j Y<2'>(\v(x)\P + \\(A j -xW \\S7vV). 

Integrating over S7j and then using the standard inverse inequality for the finite ele- 
ment function v yields 

K^n^i < c [ \v( x )\ p dx. 

By summing the above over all the nodal points Aj we obtain 

JV 



J2\v(A 3 W\{l 3 \<C [ \v\Ux, 



where we have used the fact that Qj overlaps with only a fixed number of other 
macro-elements . 

3. Maximum Principles for PI Conforming Approximations. The goal 
of this section is to establish a maximum principle for PI conforming finite clement 
approximations Uh arising from the formula (|2.3[) . This shall be accomplished by using 
a technique known as the De Giorgi's iterative method ([5]) originally developed for 
second order elliptic equations associated with maximum principles. In its essence, 
the De Giorgi's iterative technique is to estimate the set 

G(k) := {x: x£ Q,u(x) > fc} 

by showing that the measure of the set G(fc) is zero for some values of fc. The center 
piece of the De Giorgi's iterative method is the following technical lemma which can 
be proved through an iterative argument, and hence the name of the method. 

Lemma 3.1. Let <p(t) be a non-negative monotone function on [fco,+oo). 

Assume that <p is non-increasing and satisfies 

(3.1) 0( s )<f__j [<£(&)] i Vs>fc>fc , 

where a > 0, (3 > 1 are two fixed parameters and M > is a constant. Then, there 
exists a number p such that 

0(fco + p) = 0. 
Moreover, one has the following estimate 

P >M[4>{k ))^- x y a 2Pi^- i \ 



A proof of Lemma 13.11 can be found in [50] • Readers can also find more applica- 
tions of this lemma in the study of partial differential equations. For completeness, 
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we outline a proof of Lemma 13. II as follows. Let p be a real number to be determined 
later, and set 

k T = k Q + p-^:, r = 0, 1,2,---. 
It then follows from (|3.1I) that the following recursive formula holds true 

(3.2) 0(fc T+ i)< [<Kkr)] P , r = 0,1,2,.-. 



We claim that (|3.2|l implies the following 

(3.3) 0(fc T )<^M, r = 0,1,2,.-. 



with some real number r > 1 to be chosen. In fact, (|3.3p can be proved by a mathe- 
matical induction. The formula (|3.3p is clearly true with any real number r > 1 when 
r = 0. Assume that (13.31) is valid for r. Now using (|3.2[> one obtains 

IfQ0(T+l)Q 

HK+i) < — - — [<KMF 



0-1 



Now if we choose r = 1 a K& then 



From this, we see that ((33|) is also valid for r + 1 if p = M[0(fc o )P _1)/a 2^ _1) . 
Now by taking r — s> +oo in f|3 .3(1 . we see that the left limit of 4> at ko + p must be 
zero. This, together with the given monotonicity of <fr, completes a proof for the De 
Giorgi Lemma. 

Let p > 2 be any real number such that 



(3.4) p < 



-oo, d = 2, 



2d 



d>2. 



Next, we introduce a number fc* defined as follows 
(3.5) K -- 



sup max.{Ihg(x), 0}, if c > 0, 
sup Ihg(x), if c = 0. 



Assumption 1. Let the form Q(w;u,v) be given by J_?.^| ), and Uh be the finite 
element approximation of u arising from \2. 3\) . For any real number k > fc*, assume 
the following holds true: 

(3.6) Q(u h ;(u h -k)-,(u h -k)+)>0. 
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We are now in a position to derive a maximum principle for PI conforming finite 
element approximations. 

Theorem 3.2. Let Uh 6 Sh be the Pl-conforming finite element approximation 
°f HIW an d (HEj) arising from the formula 12. 3\) . Denote by I^g the interpolation of 
the Dirichlet boundary data A1.6\) that was used in the finite element formula H2.3\) . 
Let p and r be real numbers satisfying \3.J$ and 1 < r < p — 1. Assume that f G 

pr 

X,(p-i)(t— i) (Q) and the Assumption [fj holds true. Also assume that 
(3.7) c(x,u h )-^V -b{x,u h ,Vu h )>Q, VieO. 

Then, there exists a constant C — (7(0) such that 
(3-8) S u V u h {x)<K + C\\f\\ ,r 

where is given by \3.5\) . Moreover, the dependence of C = C(f2) is given by 



Proof. Let A; > A:* be any real number. Denote by tp — (u^ — k) + the positive 
part of Uh — k at nodal points. Since k > fc* and fc* is no smaller than the maximum 
value of the finite clement solution it/, on dfl, then i/j must vanish on the boundary of 
ft; i.e., 

(3.9) tp(x) e S° h . 

Thus, tp is eligible as a test function in the finite element formulation (|2 . 3[) . By taking 
v = ip in (|2.3|) . we obtain from (|2.5|) and the Assumption [T] that 

F((p) = Q(u h ;u h ,p) 

= Q(u h ;u h , {u h - k) + ) 

> Q.(u h ; (u h - k) + , (u h - k) + ) +£l(u h ; (u h - k)-, (u h - k) + ) 

(3.10) > Q(u h ; (u h - k)+, (u h - k)+). 

Using the notation p = (uh — k) + in ()3.10j) we obtain 

(3.11) {aVp, Vp) + (b • Vtp, p) + (ctp, p) = QK; tp, tp) < F(tp). 
Since the usual integration by parts implies 

(b • Vtp, tp) = -(tp, b-Wtp)- {tp, (V • b)tp), 

then we have 

{b-Vtp,tp) = -^((V-b)tp,tp). 
Substituting the above into (I3.11[) yields, 

{aVtp, Vtp) + He - iv • b)tp, tpj < F(tp), 
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which, along with the condition (|3.7p . leads to 

(aVy>, V</>) < F(tf). 

Now let G(k) be the subset of f2 where ip > 0; i.e., 

G(fc) = {T : T e T h , ¥> > for some x e T}. 

Denote by |G(fc)| the Lebesgue measure of the set G(fc). We are going to show that 
|G(fc)| = for sufficiently large values of k. To this end, we apply the ellipticity (|1.2[) 
and the usual Holder inequality to (|3.1ip to obtain 



(3.12) 



A / \Vip\ 2 dx < \\ip\\LP(a))\\f\\L<i(G(k)), 



where p > 2 satisfies (|3.4[) and q is the conjugate of p: i.e., - + - = 1. Here the L q 
norm of / was taken on the support of p for the obvious reason. Combining the usual 
Sobolev embedding with the estimate (|3.12p yields 

(3-13) imiL < C\\Vp\\ 2 L2 < C\\f\\ Lq{G{k)) \\p\\ LP . 

It follows that 

llvl|w<C||/|U. (G(fc) )<C||/|U«.|G(fc)|^, 

where r > 1 and - + - = 1 are arbitrary real numbers. The above inequality can be 
rewritten as 

Wv\\lv<c\\f\\iMk)\%. 

Now using the norm equivalence (|2.7j) we obtain 

N 

(3.14) Go^[K-fc) + (^)]P|^| <C||/||*,.|G(fc)|£. 

3=1 

It is not hard to see that G(fc) is the union of all the macro-elements ilj so that 
it/, (A,-) > k. For any p > k, one would have a corresponding set G(p). Moreover, 
if ftj C G(p), then we must have Uh(Aj ) > p > k. This implies that fl,- C G(k). 
Therefore, we have 



N 

0> - k) + (A 3 )r\tl 3 \ > Go £ [K - fcMAjOH^-, 

i =1 3=l,-,JV;i* h (^)>/!i 

>G (p-fcf £ |n 



3=1,- ,W;«h(A J )>p 



>G (p-fc) p |G(p)|. 
Substituting the above inequality into p,14[) gives 
(3.15) (p-fcf|G( P )|<G||/||^|G(fc)|^. 
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Thus, for any p > k, we have 



p 



Note that q — — and s = -Ar. Thus, 
* p— l r— l 1 

TOi<( c|l/ '';!f^ yiG(fc)i^. 

Since, by assumption, p > 2 and 1 < r < p — 1, then we have > 1. Thus, with 
4>(s) = \G(s)\, it follows from the De Giorgi's Lemma T3. II that 

(3.16) \G(d + K)\=0, 

where 

d = C2^^\vl\ e ^ l ii/n p r 

The equation (|3.16[l implies that < d + fc* on J7, which can be rewritten as 
supti/j < fc* + C2*- 1 - r |fi| 

This completes the proof. □ 



The rest of this section will establish another discrete maximum principle for the 
underlying quasi-linear second order equation when / < 0. The result can be stated 
as follows. 

Theorem 3.3. Let Uh £ Sh be the PI -conforming finite element approximation 
of (GUP and il.6\) arising from the formula $2.3\) . Let f < be any locally integrable 
function, and the ellipticity ll.ty) and the boundedness 11.3]) are satisfied. Assume 
that the Assumption [7] holds true. Then, we have 

(3.17) suptifc^) < { s ^e»^W.°). % c > 0- 
xen { s ^Pxeon hg(x), if c = 0, 

provided that the meshsize h is sufficiently small such that 

(3.18) hv<l. 



Proof. Assume that the maximum principle (|3. 171) does not hold true. We show 
that such an assumption shall lead to a contradiction. To this end, using the notation 
as given in p.5[) , we see that k* < ku = sup xgf2 Uh{x). Let fc# be the largest nodal 
value of Uh (including the nodal points on the boundary of fi) that is smaller than 
k\j- Let k be any real number such that k# < k < kja- Let ip = (v,h — k)+ G Sh be 
the positive part of Uh — k at nodal points. Since k > k# > fc* and fc* is no smaller 
than the maximum value of the finite element solution Uh on dil, then ()3.9|) holds 
true. By choosing v = f in (|2.3[) . we obtain from (I2.5[) and the assumption of / < 
that 



> F(<p) = Q(u h ; u h , <p) = Q(u h ; u h , (u h - k)+) 
(3.19) > 0.(u h ; (u h - k) + , (u h - k) + ) + 0.(u h ; (u h - (u h - k)+). 
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Now using the Assumption [T] and the notation of ip = [u^ — k) + we obtain 

£l{u h ;cp, ip) < 0, 

which leads to 

(3.20) (aVip, Vip) + (b • Slip, <p) + (ap, <p) < 0. 

Thus, we have from the ellipticity (|1.2p . the boundedness (| X . 3[) . and the condition of 
c > that 

A||V^|||2 < (aV^,V^) 
< l(b-V^)| 

(3-21) < \v\\V<p\\ L i \\ip\\m Dk ), 

where Dk is the subset of ft on which Vip ^ 0. Note that Dk is a collection of 
triangular or tetrahedral elements. It follows from the last inequality that 

(3-22) \\V<p\\ L * < i/|M| LW - 

The inequality (|3.22|) can be rewritten by using element integrals as follows 

(3.23) Yl I l V ^ 2(ff ^ v2 E / \^ dT - 

TeD k Jt TeD k Jt 

On each T C Dk, since V(f ^ 0, then ip is not a constant on T. Therefore, the 
selection of k implies that ip = at one of the vertices of T. Assume that ip{A) = 
with A being a vertex point of T . Then, we have from ip(x) = (x — A) ■ Vip that 



/ \ip\ 2 dT < h\ [ \\7ip\ 2 dT, 
Jt Jt 



where is the diameter of the element T. Substituting the above into (|3.23[) we 
obtain 

(3.24) I \^V?dT < v 2 h T I \Vtp\ 2 dT < v 2 h 2 ^ f \\7<p>\ 2 dT, 

TeD k Jt T€D k Jt TeD k Jt 

which leads to 

1 < hv. 

The above inequality is an obvious contradiction to the assumption of hv < 1 as given 
in ()3.18p . This completes the proof. □ 

4. Nodal Basis and Geometry of Finite Elements. On each triangle or 
tetrahedron T G 7~h, the finite element function v G Sh is a linear function and can 
be represented by local shape functions li — £i(x) defined as follows: (1) t L is linear 
on T, (2) £i(A(j)) = 5ij where Sij is the usual Kronecker symbol (see Fig. I4.1j) . The 
local representative property asserts that 

(4.1) v(x) =^2v(A(i))£i(x), VxeT. 

i=l 
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A(3) 




Fig. 4.1. A triangular element with acute angles 




A(2) 

Fig. 4.2. A tetrahedron with acute interior dihedral angles 

Note that the gradient of a function ip = ip( x ) 1S a vector along which the function 
ip increases the most. Thus, the gradient of the shape function ti would be parallel 
to the outward normal direction of the edge/face opposite to the vertex A(i); i.e., 

V£i = atin(i), 

where n(i) represents the outward normal direction to the edge/face opposite to the 
vertex A(i) (see Fig. 14.11 and Fig. I4.2|) . Denote by ||£|| the f 2 -length of any vector 
(elf It follows that 

on = -\\Vli\\- 

Thus, we have 

(4.2) VA = -||V4||n(i). 

The angles of the triangle AA(1)A(2)A(3) (see Fig. 14. ip can be characterized by 
using the outward normal directions n(i). For example, the angle «23 is related to 
the angle of the two normal vectors n(2) and n(3) as follows: 

a 23 = 7r-Z(n(2),n(3)), 

where Z(n(2),n(3)) stands for the angle between n(2) and n(3). Likewise, for the 
tetrahedron T as depicted in Fig. I4.2( the interior angle between the two planes 
P(A(1),A(2), A(4)) and P(A(2), A(3), A(4)) can be defined as 



= 7r-Z(n(l),n(3)). 



13 



The angle 6 is known as an interior dihedral angle. The definition of other five 
interior dihedral angles for T can be defined similarly. For simplicity, we introduce 
the following notation: 

(4.3) oy:=7r-Z(n(i),ii(7)). 
It follows from that 

(4.4) aij =w-Z(V£uVej). 

The triangle T is called non-obtuse if all the angles satisfy ay < n/2. It is said to be 
acute if < tt/2. Likewise, a tetrahedron T is called acute if each of its six interior 
dihedral angles is less than 7r/2 in radian; T is said to be non-obtuse if all six interior 
dihedral angles are no more than tt/2 in radian. For the purpose of the maximum 
principles for finite element approximations, we introduce the following concept. 

Definition 4.1. The finite element partition 7~h is called 0(h a )-acute if there 
exists a parameter 7 > such that for each element T G Tk we have onj < ■§ — 7/1", 
where a > and h is the meshsize of 7~h ■ 

5. Verification of the Key Assumption for DMP. Recall that the validity 
of DMPs as shown in Theorems 13 . 2 1 and 13 . 31 is based on the Assumption [T] which states 
that 

(5.1) Q(u h ;(u h -k)-,(u h -k) + ) >0 

for all k > fc* . The goal of this section is to verify the above assumption under certain 
conditions for the finite element partition 7ft. 

5.1. An Element-Based Approach. By an element-wise approach, we mean 
a representation of the form 0(uft; (v,h — (lift — fc)+) as integrals over each element 
T 6 7ft. To verify the assumption f|5.1[) . we shall explore conditions that make each 
element integral be non- negative. To this end, on each element T S 7ft, we use the 
local shape functions £j to represent both (uh — k)- and (tt/, — k) + as follows 

d+l 

(u h - k)-{x) = £(u„(A(»)) - fcMi(z), 
»=i 

d+l 

K - k)+(x) = ^2(u h (A(j)) - k) + l 3 {x). 
3=1 

Denote by <p — (u^ — k) + and %j) = (itft — It follows that 

0.(uh] {uh - k)-, (uh - k)+) = Q(uh] ip, tp) 

= {(aV^,\7ip)T + (h-V^,ip) T + (c^,(p)T} 
Ten 

On each element T, we have 

(aVtp, V<p)t + (b • Vip, <p)t + (cip, <f)r 

d+l „ 

= {u h {A(i)) - fc)_ - *)+ / {aV4 • + b • + c£ t £j}dx. 
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Using the angle relation (14.4)) we obtain 



Vl t ■ V£j = \\V£ 3 \\ cos{Z(V£ u Wj)) 

= ||V^||co8(7r-Oij) 
= -||V4|| ||V^||c08(ay). 



Thus, it follows from the boundedness (|1.3|) that 



/ {aVii ■ Wj + b • {Wti)lj + ctdjjdx 



Jt 



I {a||V4|| \\VtjW cos(a 2i ) - b • (Wi)£j - ci^dx 



Jt 



> [ {a\\Vli\\ ||V^||cos(a i3 -) - ||b|| ||V4|| - \c\}dx 



Jt 



> / {a||VA|| [|V^[[cos(ai 3 -)-Ai/ |[V£f|| — A^} das. 



Jt 



Assume that the element T is non-obtuse (i.e., < ctij < 7r/2). Then we have from 
the above inequality and the ellipticity (|1.2p that 



Next, we see from Taylor expansion, for onj G {po,ir/2] with po > being a fixed 
angle, there is a constant 7* > such that 



Observe that both and ||V^-|| are of size Olhj, 1 ) where hx is the size of T. 

Thus, with \T\ being the measure of T, we have 



for some A* > when the size of T is sufficiently small and ir/2 — otij > 7/1 for a large, 
but fixed constant 7. In the case of b = 0, the angle requirement can be weakened to 
7r/2 — aij > 7/1 2 . The very same argument holds true if Uh is replaced by any finite 
element function v G Sh- The result can be summarized into a lemma as follows. 

Lemma 5.1. Let v 6 Sh be any finite element function and k any real number. 
Assume that the ellipticity \1.2)) and the boundedness \1.3\) hold true. Assume also 
that the partition Th is 0(h a )-acute. Then, the following results hold true: 






>A*||V4|| IIV^'H \T 
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(i) For general b and c > 0, with a = 1, we have 

(5.2) Q(v;(v-k)-,(v-k)+)> 

a* J2 Ei(^))- fc )-i MA(j))-k)+\ iiv^h iiv^h \t\, 

Ten ift 

provided that the meshsize h for the partition Th is sufficiently small. Here 
A* is a positive number smaller than A and \T\ stands for the area or volume 
of the element T. 
(ii) For the case b = and c > 0, with a — 2, we have 

(5.3) &(v;(v-h)-,(v-k) + ) > 

A* MA(j))-k)+\ ||V4|| ||V^|| \T\, 

Ten ift 

provided that h is sufficiently small. 
(Hi) For the case of b = and c = 0, we have 

(5.4) Q(v;(v-k)-,(v-k)+) > 

A J2 £lM^«)-*)-l IK^0'))-*)+l IIViiH IIV^]] co S ( aij )\T\, 
Ten ift 

as long as each T G Th is non-obtuse. 
In other words, the Assumption^ is satisfied if the finite element partition Th satisfies 
certain angle conditions. 



B 




A 



Fig. 5.1. An interior edge shared by two elements T\ and T2. 

5.2. An Edge-Based Approach. By an edge-wise approach, we mean a rep- 
resentation of the form Q(uh] (wft — k)-, (uh ~ k)+) as integrals over macro-elements 
that share a common edge. To verify the assumption (I5.1[) . we shall explore conditions 
that make each integral on macro-elements be non-negative. To this end, we use the 
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notation if = (u^ — k) + and ip = (u;, — to arrive at 
(«/> - &)-, {uh - k) + ) = £l(u h ;ip,tp) 

d+l 



"eT h i,j=i Jt 

2 - 

_„n „_1 Jr, J 



e m „G£° 



where 8 ° denotes the set of all interior edges, A m and A are two end points of the 
edge e mn , T\ and T2 share e mn as a common edge. In Fig. I5.ll one may identify A m 

(s) 

with A, and A n with £>. Here t m is the shape function on the element T s associated 
with the vertex point A m . Thus, the validity of various DMPs can be derived if the 
following holds true 

2 . 

(5.5) £ / { aW ™ ■ Vi n ] + b • (V'&W + dx < 0. 

s=l 

In the case of Poisson problem, one has a = 1, b = 0, and c = 0. Thus, it suffices to 
have 

2 r 

(5.6) V / W« • vei^dx < 0. 

s =i Jt. 

It was known that (see for example [5]) 

V£ 1} • V6} ] d X : 



cot (a) 



and 



It follows that 



cot(/3) 



E / • w « w = 



cot(a) cot(/3) 
2 2 

sin(a + j9) 
2 sin a sin /3 ' 



and (|5.6[) holds true if and only if a + /3 < tt. 

A similar, but more complicated, analysis can be conducted for tetrahedral ele- 
ments; this is left to readers with interest and curiosity on DMPs for Poisson problems 
in 3D. 
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